Intermolecular Effect in Molecular Electronics 
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We investigate the effects of lateral interactions on the conductance of two molecules connected 
in parallel to semi-infinite leads. The method we use combines a Green function approach to 
quantum transport with density functional theory for the electronic properties. The system, modeled 
after a self-assembled monolayer, consists of benzylmercaptane molecules sandwiched between gold 
electrodes. We find that the conductance increases when intermolecular interaction comes into play. 
The source of this increase is the indirect interaction through the gold substrate rather than direct 
molecule-molecule interaction. A striking resonance is produced only 0.3 eV above the Fermi energy. 

PACS numbers: 72.80.Le, 85.65. +h, 73.40.Cg 
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Modern molecular electronics began in 1974 when Avi- 
ram and Ratner 1] proposed a rectifier based on an asym- 
metric molecular tunneling junction. The size of molecu- 
lar electronic devices - on the nanometer scale - provides 
advantages in cost and efficiency, while synthetic modi- 
fication of molecules may enable manipulation of their 
electrical properties. All these advantages promote the 
possible use of such systems in future molecular electronic 
technologies. In fact, several fundamental devices have 
recently been demonstrated 0,111 DH3 In many of these, 
self-assembled monolayer (SAM) systems were used as 
they are a controllable step toward single molecule elec- 
tronics. 

Molecules thiolated on a gold surface comprise a fre- 
quently studied SAM. For instance, scanning tunneling 
measurements show that conjugated molecules with one 
methylene group between the sulfur and aromatic ring 
have a commensurate (v^x y/3)R30° structure Q- Sul- 
fur is chosen because it has two valence electrons, allow- 
ing bonding to both molecule and surface. Among other 
possibilities, oxygen has high electronegativity leading to 
a high dipole barrier, while Se and Te result in lower con- 
ductance as suggested by calculations on a simple phenyl- 
dithiolate molecule [7J. At high coverage, the SAM forms 
via the strong S-Au bond 0, EJ; a strong chemical 
bond is more likely to generate stable, reproducible junc- 
tions and avoid a large charge-transfer induced barrier 
which may mask the electronic signature of the molecule. 

The conductance of a SAM may involve both intra- 
molecular and inter-molecular transport. In the case of 
densely packed conjugated molecules 0, 0, 0, the 
lateral interaction through delocalized 7r-orbitals could 
mask the individual contribution of a presumably iso- 
lated molecule. Thus, when two molecules with con- 
ductance G\ and Gi are sandwiched in parallel between 
two leads, the resulting conductance is not simply 
equal to G\ + G%. Such considerations cast doubt on the 
common procedure whereby the conductance of a single 
molecule is deduced from Gtotai by simply dividing by the 
number of molecules attached to the leads. Conversely, 
theoretical calculations on a single molecule may not be 
relevant for comparison to SAM experiments. 

For parallel molecular wires that connect directly to 



the electrodes, indirect interaction was found by Yali- 
raki and Ratner |14|. They set the interchain hopping 
to zero, and saw an increase in conductance, in the ab- 
sence of any coupling between two molecules. A linear 
superposition law has been found using a tight-binding 
description |l5j . The intensity of coupling as a function 
of equilibrium site distance was taken as evidence for in- 
teraction through the Au surface. Lang and co-workers 
[T(| calculated the low bias conductance of cumulene on 
jcllium as a function of molecule-molecule distance. The 
range of intermolecular distance considered, however, is 
3.5-6.5 a.u. - less than 9.4 a.u., the distance between the 
closest adsorption sites of Au (111) surface. 

Since aromatic molecules are widely used in molec- 
ular electronics, we study the effect of intermolecular 
coupling on the conductance of two benzylmercaptanes 
(BM) sandwiched between semi-infinite gold leads. The 
geometry of the molecules (Fig-QJ is that of possible crys- 
talline structures for the infinite SAM. The suggested 
herringbone packing of a BM-SAM [6j has two kinds 
of mutual structures: displaced tt stacking and verti- 
cal packing. Here we focus on the displaced tt stack- 
ing because of the shorter distance between delocalized 
7r-bonds. Our density functional theory plus Green func- 
tion approach has been described elsewhere . We used 
the Kohn-Sham method in the generalized gradient ap- 
proximation for the electronic structure calculation, with 
numerical atomic orbitals (double zeta plus polarization) 
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FIG. 1: (Color online) Structure of (a),(b) single benzylmer- 
captane molecules at each of two bridge sites on a Au(lll) 
lead and (c) the corresponding double molecule system. Each 
Au layer contains 9 atoms. Brown, white, blue and yellow 
denote Au, H, C, and S atoms, respectively. 
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TABLE I: The initial and relaxed site, the adsorption energy 
with respect to the fee-bridge case, the height of the sulfur 
atom from the gold surface, the average of the three nearest 
neighbor Au-S distances, and the S-C bond length and angle 
9 to the surface normal. (81 fc-points are used, yielding a 
converged adsorption energy |2(i|.) 
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as the basis set, by means of the SIESTA program 18]. 

Fig. Q] shows the structure of one or two BM moleclues 
adsorbed on the gold surface at bridge sites. The dis- 
tance between the two equilibrium adsorption sites for 
the BM moelcules is 9.4 a.u. Unlike in a break junction 
[lfij . each molecule has only one S-Au contact. A second 
contact for conductance measurements is provided by a 
STM tip contacting the top hydrogen atom. We start by 
considering single adsorbates in high symmetry sites on 
the Au(lll) surface - top, bridge, fee hollow, and hep 
hollow. According to our calculations, the hep and top 
sites are not local minima; this agrees with methylthi- 
olate results [2(j showing that the hep site is unstable 
at high coverage. As shown in Table |U the tilt angle of 
the relaxed fee conformation is 7° less than that of the 
relaxed bridge site. This implies more steric repulsion 
between the surface and the C-S back bond. In the end, 
the bridge site is favored by 0.16 eV. 

The conductance and charge transfer changes due to 
intermolecular coupling are shown in Table ITT1 The small 
difference between Q\ and Q2, or G\ and G2, is due to 
the different position on the surface. We see that the 
change in charge transfer due to the insertion of a second 
molecule is very small. In contrast, the conductance in- 
crease is large, 65% and 58% for adsorption at the bridge 
and fee sites, respectively. This demonstrates clearly that 
the conductance of a N-molecule SAM is not simply N 
times the conductance of a single molecule. 

The charge transfer is essentially a local issue because 
the sulfur atom does not connect directly to the conju- 
gated aromatic ring. Therefore, the charge transferred to 
the sulfur cannot be efficiently shared with the core part 
of the molecule (aromatic ring). The methylene group 



TABLE II: Charge transfer from the leads to molecule 
and conductance of single and double molecules adsorbed at 
bridge or fee sites. The subscript "1" denotes a single molecule 
at the central adsorption site, "2" a single molecule closer to 
the edge of the lead, and "1,2" two molecules. AQ is de- 
fined as (Qi,2 — Qi — Q2)/{Q\ + Q2) and AG similarly. The 
conductances are in units of 10~ 4 x 2e 2 /h. 



Site 


Qi 


Q2 


Qi,a 


AQ(%) Gi 


G 2 


G1.2 


AG(%) 


bridge 


0.202 


0.071 


0.266 


-5.5 2.6 


1.7 


7.1 


+65 


fee 


0.199 


0.110 


0.292 


-2.5 3.5 


7.1 


16.8 


+58 



between the sulfur and the aromatic ring does provide a 
support tool so that the BM will be straight up when ad- 
sorbed; on the other hand, it destroys the conjugate char- 
acter of the molecule and so blocks the charge transfer. 
This trade-off between mechanical and electrical proper- 
ties is a fundamental dilemma for molecular electronics: 
the inserted methylene group makes the molecule rela- 
tively insensitive to environmental changes - crucial for 
reproducibility - but its conductance is low. 

We now use the equilibrium transmission, T(E), and 
the density of states projected on the BM molecule(s), 
PDOS, to further analyze our system; see Fig. [21 The 
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FIG. 2: Transmission (dashed line) and PDOS (solid line) 
of BM molecuels adsorbed at bridge sites: (a) single molecule 
at site 1, (b) single at site 2, and (c) double molecule. The 
inset in (c) shows surfaces of constant local density of states 
(LDOS) at the HOMO (left) and LUMO (right) energies. 
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FIG. 3: (Color online) (a) A surface of constant LDOS at 
0.3 eV for the double molecule superposed on a ball-and-stick 
model. Note the involvement of 7r-orbitals on the central (site 
1) BM. (b) Difference between LDOS of double molecule and 
that of single site-1 case. The blue and red denote positive 
and negative surface values, respectively. Brown, white, blue 
and yellow denote Au, H, C, and S atoms, (bridge site) 
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FIG. 4: Difference of DOS between double molecule case and 
sum of both single molecule cases projected on the aromatic 
rings (dashed) and the S-Au interface (solid). For the latter, 
the S atom plus the two nearest layers of Au are used. The 
change in PDOS near the Fermi energy clearly comes from 
the indirect interaction through the Au. (bridge site) 



coupling to the continuum of metal states leads, of 
course, to shif ting and broadening of the molecular en- 
ergy levels 0,123- Note that the curves for the two sin- 
gle molecule cases are very similar [panels (a) and (b)]; 
the difference comes only from different location with re- 
spect to the edges of the lead. The HOMO-LUMO gap 
is 3.8 eV in an isolated BM molecule calculation, and 
the same energy scale remains when connected to leads. 
From the inset of Fig. Etc), we see that the HOMO is 
localized at S atoms, which contribute low transmission. 
On the other hand, the LUMO is a delocalized 7r-orbital 
and leads to large T. In the gap, the transport process 
is non-resonant tunneling yielding low transmission. 

Both T(E) and the PDOS change upon going to a two 
molecule system (Fig. |2J). A very striking feature ap- 
pears: a narrow resonance about 0.3 eV above the Fermi 




FIG. 5: Surfaces of constant LDOS at the Fermi energy su- 
perposed on models of the molecules, (a) Single BM molecule 
at site 1. (b) Single molecule at site 2. (c) The Au lead with- 
out BM molecules, (d) Double molecule, (e) Difference be- 
tween single molecule at site 1 and bare lead, (f) Difference 
between double molecule and bare lead. Note the electron 
deficient areas around the S-Au bonds; with two molecules 
present, this deficiency is decreased thus leading to better 
transport properties. Blue and red denote positive and neg- 
ative values, respectively. Brown, white, blue and yellow de- 
note Au, H, C, and S atoms, respectively. The energy window 
is ±0.02 eV around the Fermi level, (bridge site) 



energy. This feature appears in both the fee and bridge 
site cases, but is not present for an isolated double BM 
system. Such a feature is relevant to molecular electron- 
ics since in the I-V curve at low bias it would cause a 
region of negative differential conductance. 

To show the nature of the resonance, we plot in Fig. 
Fig. 01a) a surface of constant local density of states 
(LDOS) for the energy 0.3 eV. The resonance is clearly 
made from LUMO related anti-bonding 7r-orbitals on the 
central BM molecule. The difference between the LDOS 
for the double molecule and that for the single site-1 case 
is shown in panel (b). An enhanced electron density ap- 
pears at the sulfur connection when molecule 2 is added. 
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The better S-Au coupling that this causes is apparently 
responsible for the resonance peak at 0.3 eV. 

Our results show clearly the large changes that inter- 
action among the molecules of a SAM can induce. Yet, 
the question remains as to what brings about those al- 
terations - direct molecule-molecule interaction between, 
for instance, the aromatic rings, or indirect interaction 
through the gold surface? 

We answer this question by looking at how the DOS 
changes from the single molecule to double molecule 
cases. Fig. 0] shows the change of the DOS projected 
on either the aromatic rings or the S-Au interface. Large 
fluctuations take place at the S-Au interface but not on 
the rings. The flucutations occur both at energies near 
the edges of the HOMO-LUMO gap, which is natural 
from shifting of levels due to changed coupling, and near 
the Fermi energy. Thus, the dominant interaction be- 
tween the molecules is indirectly through the Au leads, 
and this controls the transport behavior 23]. 

To further understand the changes in conductance at 
the Fermi energy, we focus on the LDOS at that energy 
(see Fig.|SJ). It is known, of course, that sulfur bonds well 
with nearby gold atoms, but we see only a small LDOS 
around the Fermi energy [panels (a)-(c)]. Thus, this type 
of bonding does not provide a good conduction channel. 
This is an example of a structurally favorable anchoring 
group which provides reproducible measurements but is 
not a good choice for a transparent contact. The bottle- 
neck of molecular electronic devices is, in fact, this con- 
nection between leads and molecule. It depends largely 
on the choice of lead material and head group atoms. 
For better performance, the coupling should not only be 



strong far from Fermi energy. In our case the LDOS does 
not spread to the core molecule from the lead, because 
of either the large gap of the molecule or the methylene 
group blocking the conjugation of the system. 

Fig. 0e) tells us that when there is only one BM 
molecule, a density deficient region forms around the S- 
Au contact at the Fermi level, another demonstration 
of the poor conducting properties of this bond. On the 
other hand, Fig. Of) shows that when two molecules are 
present, there is an increase of LDOS in between them. 
The S atom tends to push electrons at the Fermi energy 
away from it; however, when two molecules get together, 
the mutual pushing in some sense eliminates the electron 
deficient region caused by a single S atom, and so yields 
a better connection with increased conductance. 

Our calculations thus provide a concrete example of 
weakly interacting parallel molecules adsorbed on a gold 
surface, one of the scenarios investigated by Yaliraki 
and Ratner [14| with a model Hamiltonian. The direct 
molecule-molecule interaction can be neglected because 
of the large HOMO-LUMO gap and limiting states at 
the Fermi level. It is the indirect interaction through 
the gold substrate that modifies the conductance. Close 
packed molecules can change the electron density at the 
Fermi energy of gold leads near the surface, and this ef- 
fect extends to more than four layers of gold atoms before 
decaying. Thus, a nearby molecule could have a dra- 
matic effect on the conduction due to a modification of 
the LDOS of gold leads. 

This work was supported in part by the NSF (DMR- 
0103003). 
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